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ABSTRACT 

Strong lensing provides popular techniques to investigate the mass distribution 
of intermediate redshift galaxies, testing galaxy evolution and formation scenarios. It 
especially probes the background cosmic expansion, hence constraining cosmological 
parameters. The measurement of Einstein radii and central velocity dispersions in¬ 
deed allows to trace the ratio Dg/Dis between the distance Ds from the observer to 
the source and the distance Dis from the lens to the source. We present an improved 
method to explicitly include the two-component structure in the galaxy lens modeling, 
in order to analyze the role played by the redshift and the model dependence on a nui¬ 
sance parameter, (Fe, which is usually marginalized in the cosmological applications. 
We show how to deal with these problems and carry on a Fisher matrix analysis to 
infer the accuracy on cosmological parameters achieved by this method. 
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1 INTRODUCTION 

Dozen years after its unexpected and somewhat serendip¬ 
itous discoverty, the accelerated expansion of the universe 
is taken for granted due to the flood of dat a from different 
astrop hysical probes confirming it (see, e.g.. lWeinberg et al.l 
ll2012ll for a not updated yet exhaustive review). Although 
the spatially flat concordance ACDM model, made out of 
a cosmological constant accounting for ~ 70% of the en¬ 
ergy budget and responsible of the cosmic speed up, is in 
full agreement w ith observations (as summarized, e.g., in 
lAde et ^ ll2015h . where the latest datasets are considered), 
it is f ar from free of any conceptual and theoretical prob¬ 
lems (ICarrollll200ll : [Pemianski fc Piedipalumbcill2015l l. This 
motivated the search for alternative models generally re¬ 
ferred to as dark energy (DE) ones and characterized by 
the presence of a scalar field with a negatiye equation of 
state (EoS) denoted as w = p/p, with (p, p) the pressure 
and energy density of this leading component. While a huge 
number of papers has addressed the problem of what DE is 
with proposal running fro m self interacting fields to modi fied 
gravity theories (see, e.g.. lAmendola fc Tsuiikawal (l2015ll for 
a recent textbook), observations are mainly aimed at con- 
straining the DE EoS. Under the simple yet efficient CPL 
dChevallier fc Polarski|[^0ll : lLinderll2003li parameterization, 
it is considered w = wo + Wa{l — a), with a = 1/{1 -\- z) the 
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scale factor and z the redshift, so that the Holy Grail of ob¬ 
servational cosmology has nowadays become to narrow down 
as much as possible the range for the {wQ,Wa) parameters. 
As a consequence, the efficacy of an observational probe is 
presently quantified in terms of the Figure of Merit (FoM) 
defined as the inverse of the area delimiting the 95% co nfi- 
dence range in the (wo,Wa) space (I Albrecht et al.|[2006l i. 

Gravitational lensing has soon been pointed at as one 
of the most promising tools to investigate the nature of the 
mysterious DE. In its journy from the source to the observer, 
the light is deflected by both a dominant mass concentra¬ 
tion (a galaxy or a cluster acting as lens) along the line 
of sight and the cumulative effect of the large scale struc¬ 
ture. As a consequence, gravitational lensing probe both the 
background expansion and the growth of structures. Gos- 
mic shear tomography, that is to say the power spectrum of 
shear as m easured from galaxies separated in different red¬ 
shift bins (iBartelmann fc Schneide^l200ll : iHoekstra fc Jairj 
l2008ll . has emerged as one of the most promising tools to 
constrain the DE EoS, so motivating the interest in on¬ 
going (DES, KIDS) and future (Euclid, LSST) dedicated 
surveys. While shear tomography is based on the weak 
regime of gravitational lensing (where lensed galaxies are 
subject to tiny modifications detectable only on statisti¬ 
cal grounds), also statistics based on strong lensing fea¬ 
tures is not less important. In this regime, the light emit¬ 
ted by the source passes close to a massive lens causing 
the formation of multiple images or spectacular Einstein 



























2 V.F. Cardone et al. 


rings or radial and/or tangential arcs. The distribution of 


angular senaration between double images dLi & Ostriker 

2002 

ilLooes Kr. Millerll2004l:lMitchell et alJ 

2005; 

Zhang et al. 

200 c 

1, arc statistics dBartelmann et aid 

to 

o 

o 

Dalai et al. 

2004 

: iMeneehetti et al. 2013: Boldrin et al 

2015ll and 

Einstein rings nronerties (iBroadhurst & Ba 

rkanal l2008l: 

Zitrin et alJl2012l: IWaizmann et alJ 20141) 

have been used as 


tools to constrain cosmological parameters and the astro- 
physical properties of the lens population. 


Galaxy scale strong lenses have been also proposed 
as a tool to probe the background expansion. The basic 
idea is usually to infer the distance ratio Dis/Ds between 
the angular diameter distance from the lens to the source 
and from the observer to the source, and then use a suf¬ 
ficiently large sample to trace how this quantity evolves 
with z. To this end, one typically relies on measurements 
of both the Einstein radius Re and the lens velocity dis¬ 
persion (tq. Under the assumption of isothermal sphere, one 
gets Re = 4:'K{Dis/Ds){asis/c)'^, so that a measurement 
of (RejO-q) immediately gives the distance ratio Dis/Ds 
provided one sets crsis = . ferro, being fs a nuisance pa¬ 
rame t er to marginalize over (Biesiada et al.ll2010l : ICao et al] 
I2OI2I : [Piorkowska et al.ll201S 1. Although quite simple, such 
a method however rests on some assumptions which, even if 
reasonably motivated, should nonetheless be verified. First, 
it is assumed that all the lenses can be modelled as singu¬ 
lar isothermal sphere. Second, the nuisance parameter /_b is 
taken as a sort of universal constant, while its actual value 
can change from one system to another depending on both 
the lens properties (redshift, mass, size) and how (tq is mea¬ 
sured (weighted within a circular aperture of fixed physical 
radius or through longslit spectroscopy). Here we follow a 
similar approach, but weakening the strong assumptions of 
the standard approach in order to avoid biasing cosmological 
parameters because of incorrect lens modelling. 

The plan of the paper is as follows. In Sect. 2, a step - by - 
step derivation of the method is presented underlying which 
assumptions are made and how we deal with them. Sect. 3 
is devoted to validate the proposed technique showing that 
its basic assumptions are well motivated. We then present, 
in Sect. 4, a Fisher matrix analysis to forecast the errors 
on cosmological parameters for different configurations of 
future datasets. A summary of the results and a discussion 
of future perspectives are given in the concluding Sect. 5. 


with (M*,M„ir) the total stellar and virial masses, respec¬ 
tively, and f^{R) is a dimensionless function depending on 
the density profile. For the stellar component this reads 


t^AR) 


r(2n)-r[2n,6„(i?/i?e//)^/"l 

r(2n) 


( 2 ) 


with n the slope of the Sersic model, bn a constant set so 
that half of the tota l luminosity is contained within the ef¬ 
fective radius Reff (|Ciottilll99ll l. We leave for the moment 
the fj.DM (R) function undetermined in order to be as general 
as possible. 

Provided that the lens and source redshifts {zi,Zs) are 
known, one can get an estimate of the projected mass by 
simply measuring the Einstein radius and then using 


M°noj{RE) = ETscritizi, Zs)Re = 


AR%{arc) DiDs 
4G(206265)2 As 


(3) 


where we have used the definition of critical surface density 
^crit Al , ^s) ■ 
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crit — 


Ds 

47rG DiDis ’ 
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being c the speed of light and {Di, Da, Dis) the angular di¬ 
ameter distances to the lens, the source and between lens 
and source. Note that, in Eq.([3ll, f?£:(nrc) = Re{‘ 206265/Di) 
is the Einstein radius in arcsec since this is the quantity di¬ 
rectly measured from the data in a model independent way. 

The second quantity of interest is the aperture velocity 
dispersion, i.e. the velocity dispersion luminosity weighted 
within a circular aperture of radi us flap. Approximat i ng the 
deprojected Sersic profile with the lPrugniel fc SimierJ (Il997fl 
model, it turns out that 


_2 GM* bA^-^’'^ri2n) 

Raff ne>>«r[n{3-pn)] ^ 

'L-k {^Rap 5 -j- j Adk (^Rap Ref f j P £) jVf ) 

r(2n) — r[2n, bniRap/ReffY^'^] 

where pn is given in iMarauez et al.l (1200 il l, pdm collec- 
tively denotes the ha lo model parameters, and we refer to 
ICardone et al.l ll2009ll and ICardone &: Tortoral (I2OIOI I for the 
details of the dimensionless T{x, p) functions. For dimen¬ 
sional reasons it is convenient to define an aperture mass 
as 


2 DISTANCE RATIO 

Galaxy scale strong lenses are typically used to constrain the 
dark matter halo and its interplay with the properties of the 
stellar component. To this end, two observational probes are 
used, namely the projected mass Mproj within the Einstein 
radius Re and the central velocity dispersion (Tq. Assuming 
the galaxy may be modelled as the sum of a stellar part 
modelled with a Sersic (1968) profile and a dark matter halo, 
we gelQ 

MlAjiRs) = Mi,fii,{RE) 


1-b 






-Pdm{Re)p*{Re) , (1) 


^ Unless otherwise stated, length quantities such as the Einstein 
radius are meant to be in linear units. When this is not the case, 
we will explicitly give the angular units adopted. 
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= M* 
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G ne^"r[n(3 — Pn)] 

T* {Rap , n) -(- (Afnir / Af* M {Rap , tl, Raf f , POM ) 

r(2n) -r[2n,bn{Rap/RaffA^A 


( 6 ) 


which can be straightforwardly estimated from measurable 
quantities as 


t^-oba _ .Re//(fflrc) DlfJap 

206265 G ’ ^ ’ 

with Raff {arc) the effective radius in arcsec obtained from 
the fit to the surface brightness profile. 

It is now only a matter of algebra to first solve Map’’ = 
Map with respect to M* and then insert the result into 
Mpnaj{RE) = Mpnnj{RE) to finally get 

D{zi, Za) — Doha X De , 


( 8 ) 
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which represents the basic relation of our method relating 
the theoretical distance ratio 

'D{zi,Zs) = Ds/Dis (9) 

on the left hand side with the {T>obs, tFs) quantities on the 
right hand side defined as follows 

_ 4 X 206265 i?e//(arc) /crap\^ ne'’"r[n(3 — pn)] 

RE{arc) RE{arc) V c / 

V[2n,b„{Rar,/ReffY''^]\ 

r(2n) / 

V[n{Z-prt),hCRap/ReffY''^] 
r[n(3 - pn)] 

j. _ + {Myir/Mi,)\p,DM{RE\'PDM)/P-k{RE\n)] 

'R*lRap^^^ (-^ap 5 n j , PPJJVT ) 

Eq.((8| is formally similar to the one used in previous sim¬ 
ilar methods expressing the distance ratio in terms of a 
quantity depending only on measurable quantities (given by 
'Dobs) and an unknown nuisance parameter (referred to as 
Fe here). However, we improve on the classical approach 
in two ways. First, we have not assumed the galaxy to be 
a one component system, but rather modelled it as it ac¬ 
tually is, i.e. the sum of a stellar parl[3 and a dark matter 
halo. Second, we have now an analytical expression for the 
nuisance parameter, so that we can investigate whether it is 
the same for all lenses or rather it is a function of both the 
lens properties and the halo model. 


3 THE NUISANCE PARAMETER 

To infer the distance ratio from strong lensing data, one 
generally assumes that the /b parameter is a constant inde¬ 
pendent on the lens and marginalizes over it. This is actually 
a zero order approximation related to the properties of the 
assumed singular isothermal sphere model used to describe 
the total mass density. On the contrary, with our Fe param¬ 
eter in Eq. (HB) we account for the two components nature 
of the lens and explicitily specify the halo model. 

A careful analysis of the different terms entering its 
definition shows that Fe depends on (i) the mass ratio 
Mvir/Mi,, (ii) the length ratio R- 2 /Reff (with R -2 the ra¬ 
dius where the DM logarithmic slope equals -2), and (hi) 
the halo concentratioi^j Cvir = Rvir / R-2- Actually, the halo 
model can be more complex, so that {M^ir, Cvir) are not 
enough to fully assign it, and more parameters are needed, 
which makes Fe depend on them too. Moreover, since the 
mass and length ratios may both change with the redshift 
as it is suggested by, e.g., the redshift dependence of the 

^ By adopting a single Sersic profile, we are implicitly assuming 
that the lens is an early - type galaxy, otherwise one should split 
the stellar term in the sum of a Sersic bulge and an exponential 
disk. This is not a serious limitation since most of the lenses are 
indeed intermediate ellipticals. 

® No matter which halo model is adopted, it is always possible to 
define the concentration provided that the logarithmic slope 
is not constantly equal to -2, which is the case for the singular 
isothermal sphere only. 


Cvir - Mvir relation as well as of the Re// - Af* scaling law, it 
turns out that Fe is redshift dependent. 

These considerations strongly point against the stan¬ 
dard assumption to consider the nuisance parameter as a 
constant, which implies that one should fit for Fe separately 
for each lens, hence making the method unusable. However, 
a possible way out could be found if Fe was constant within 
a certain scatter once the lenses were binned according to 
some suitably chosen properties. To investigate whether this 
is possible or not, we proceed as schematically pointed out 
below. 

(i) We set the lens redshift zi sampling from a uni¬ 
form distribution over the range (0.1,2.1) and then set 
Zs = Czi with g randomly extracted from the range (1.1, 3.1). 

(ii) We use t he redshift dependen t galactic stellar mass 
function used in iFontana et al.1 ll2006h to set the lens stellar 
mass M* and sample the effective radius Ref / from a 
Gaussian distribution centred on Rref{Mi,,z = 0)(1 + z)°‘ 
and a 10% scatte r. The estimation of Rref{Mt,z = 0) is 
obtained following IShen et al.l (l2003li . while t he scaling with 
redshi ft has been chosen in accordance with iTruiillo et al.l 
1 I 2 OO 6 D considering the value of a for massive galaxies 
(M* > 3 X 10^° Mq). 

(iii) The halo virial mass is set s olving the vs 

Moir relation of iMoster et al.l (l2010l') with respect to Mvir, 
so that we can finally fix the concentr a tion s ampling from 
the Cvir-Mvir relation of iDuffv et al.l (120081') . taking also 
into account its scatter. 

(iv) Adopting an NEW dNavarro et al.|[l9^ . ll997l ') model 
for the DM halo, we can now compute both the Einstein ra¬ 
dius Re and the central velocity dispersion cto = (Tap{Rap = 
Reff/8). We then retain the simulated lens only if its input 
and output quantities are reasonable. To check whether this 
is the ca se, we compare them to the similar quantities for the 
SLAGS d Auger et al.ll20(3^ sample of 85 intermediate red¬ 
shift lenses with well measured values of [zi, Zs, M*, Rb, uq). 

Using the above recipe, we can generate a sample of 
« 10000 lenses, with a distribution in [zi, Zs, M-,,, Re,(to) 
similar to those of the real SLAGS, and confidently use it 
to investigate how the nuisance parameter, Fe, scales with 
the lens properties. Moreover, since we are interested more 
in the distance ratio estimate than in Fe, we wonder what 
could be the bias on D{zi,Zs) caused by using the mean 
value Fe ~ (Fe), averaged over a subsample of 1000 lenses 
selected according to a given criterion, instead of Fe- To 
this end, we define the distance ratio estimator 

■D = Dobs X {Fe) (12) 

and consider the quantity 

A = (R - V)/D = {Fe - {Fe))/Fe 

as a measure of the relative error due to the distance ratio 
estimate. Since {Fe) is the average value over the lenses in 
a given bin (so that we can trace how A evolves with the 
bin centre), the question is now which parameter to use as 
binning quantity. 

As a first case, we consider the most obvious choice and 
divide the simulated sample in 50 redshift bins. It turns out 
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Figure 1. Nuisance parameter Te and scatter A as functions of redshift zi, stellar mass logM*, central velocity dispersion ctoi and 
Einstein radius Re for the 1000 lenses subsample. Results for the whole 10000 lenses sample are qualitatively the same. 


that, while {Te) increases with zi, its relative scatter only 
has a negligible variation, and that the root mean square 
value is Arms — 20% independently of zi. Moreover, we have 
verified that, for lenses in the same bin, A has a negligible 
correlation with the source redshift, the lens size and mass, 
and the halo model parameters. 

We repeat this test, binning with respect to the stel¬ 
lar mass, central velocity dispersion, and Einstein radius. 
We find that we are unable to get satisfactory results. In¬ 
deed, binning with respect to logM* increases Arms and 
introduces a significant increase with the stellar mass for 
logM* > 11.2, with values going as large as ~ 50%. Sim¬ 
ilarly large values are obtained when binning with respect 
to (To although a definite trend is not found, while binning 
in Re introduces an increasing trend for Re < 1.5 arcsec 
followed by a non monotonic variation leading to Arms — 
30 — 40% (but the low statistics at large Re values makes 
these numbers not fully reliable). 

We therefore conclude that the estimator GU is reli¬ 
able, once the lenses are divided in redshift bins. In this 
case, one can then rely on "D assuming that the nuisance pa¬ 
rameter is the same for all the lenses in the same bin within 
a percentage scatter of ~ 20%. 

To further validate this conclusion, we have repeated 
the full analysis changing the assumptions in the simulation 
pipeline. We have indeeed considere d a different Cyjr - M^ir 
relatio n either by adopting the one of iMunoz - Cuartas et al.l 
1 I 2 OIIII or by changing the halo model from NEW to Einasto 
(1965, 1969) profile. Although the Te values are different, 
the results on Arms are qualitatively and quantitatively sim¬ 
ilar. In particular, we confirm the above conclusion that 
binning in lens redshift guarantees that Arms — 20% in¬ 
dependently of the bin centre and the stellar and DM halo 
properties of the lenses. 


4 FISHER MATRIX FORECAST 

In order to investigate whether and under which conditions 
the method proposed works in constraining the cosmolog¬ 
ical parameters, we perform a Fisher matrix analysis. The 
elements of the Fisher matrix are given by 


F ■ — 


d^C 


dpidpj 


p=pfid 


being Pi the i-th parameter and Pfid the fiducial values. To 
define the likelihood function C we split the sample in A/i, 
bins and let A4 be the number of lenses in the fc - th bin. We 
can compute the likelihood for lenses in this bin as 


= 


+ 


^ ^ In (ud-b (j|) 



d=l 


2 


T^obs,d{'^E) k Pc) 

V^d + ^l 


(13) 


where pc is the set of cosmological parameters to be con¬ 
strained, while {J-E)k is the average Te value for the k — th 
bin which we marginalize over. In order to compute the like¬ 
lihood (and its derivatives) we need to assign the quantities 
entering Eq. (II3J. The 'Dobs,d may be easily computed for 
each lens using Eq. dini; the errors {ad,crE) deserve, on the 
other hand, some more words. First of all, ad is the error on 
'Dobs,d{TE)k- Since we are assuming that {RE)k is constant, 
we can then simply write ad = {TE)kaobs, where a naive 
propagation of errors is sufficient to get 


^ohsj'Rohs^d ^E ^e// S” ‘^^0 ■ 

Here {eE,Eeff,Eo) are the relative uncertainties on 
{Re , Ref f , ao) , respectively. We will make the simplifying 
assumption that these quantities are the same for all the 
lenses in a given bin and explore different values for them 
(see later). The second term aE is the systematic error in¬ 
duced by replacing the actual Re with its average value and 
can be set as as ~ Arms {Re) k- Based on the simulations 
described in the previous section, we will set Arms = 0.2, 
while {RE)k is computed by binning the lenses according to 
the above chosen strategy. 

Assuming that the lens redshift has been measured with 
a precision much smaller than the bin width, the total like¬ 
lihood is simply: 
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j^b 

CPc, Pe) = £fc(Pc, {tFE)k) , 
k=l 

where we collectively denote with p_B the set of {J-E)k values. 
As a consequence, the total Fisher matrix will be the sum 
of the Fisher matrices for each bin. It is convenient to first 
compute the Fisher matrix for the k-th bin, then marginal¬ 
ize over {J-E)k, and finally sum the marginalized matrices 
to get Fc, i.e., the total Fisher matrix for the cosmological 
parameters only. 


4.1 Input quantities 

In order to compute the Fisher matrix, some preliminary 
quantities must be set. First, we choose the fiducial cos¬ 
mological model. We consider three different spatially flat 
models setting 

r (0.306,-1.0,0.0) 

(O.Mj'Wo^Wa) = I (0.306,-0.95,0.0) 

[ (0.306,-0.90,-0.16) 

which we will refer to as ACDM, Quiessence, Thawing, re¬ 
spectively. Note that we have used the same matter density 
parameter for all three models, setting its value to the recent 
Planck estimate although one should change it according to 
the DE EoS parameters in order to get the best match to 
the data. However, we are here only interested in exploring 
realistic models so that we do not worry much about the ex¬ 
act value of the parameters. We nevertheless want to stress 
that these cases span the range of q uintessence models still 
allowed by the data (lLindejl2015al l. These parameters as¬ 
sign the dimensionless Hubble parameter E(z) = H{z)/Hq, 
reading 


E{z) = tluil + zf 

+ (1 - F!m)(1 + ^)3(i+“o+“«) exp (-^) , (14) 

which enters the distance ratio given by 


Eizi, Zs) 
with 


Ds ^ xCs) 

Dis xCs) - xCi) ’ 


XC) = [ E-\0d(; 

Jo 


(15) 


(16) 


the comoving distance. Note that Eg. (1151) only holds for spa¬ 
tially flat models. We also stress that the Hubble constant 
Ho drops out from the ratio of distances thus reducing the 
number of cosmological parameters. 

A key role in determining the accuracy of the con¬ 
straints is played by the errors on the observable quanti¬ 
ties. These enters Eq. (na through the statistical term aobs 
and the systematic one (t_b. As already said above, we set 
(je = Arms{EE)k, since this is the scatter in the distance 
ratio introduced by the approximation done when using the 
estimator instead of the correct one in Eq. ®. On the 
other hand, the statistical uncertainty aoba depends on the 
relative errors on (Re, Reff, fro). We will set e^ff = 0.015 as 
typically found fitting HST - like data on lens surface bright¬ 
ness profile, while we try three different conhgurations for 


Table 1. Marginalized constraints on (r2jvf,tuo) for ACDM, 
Quiessence and Tahwing models (top, centre and bottom part 
of the table, respectively) for different sample and error config¬ 
urations. The labels Ik and lOfc refer to the sample containing 
1000 and 10000 lenses, respectively. We remind the reader that 
we set £e// = 0.015 for all the cases. 



£0 

GlkC^M) 

o'ife(n)o) 

GlOkC^M) 

Giok (rco) 

0.05 

0.10 

0.026 

0.195 

0.008 

0.060 

0.01 

0.10 

0.024 

0.184 

0.008 

0.056 

0.01 

0.05 

0.019 

0.140 

0.006 

0.042 

0.05 

0.10 

0.027 

0.191 

0.008 

0.058 

0.01 

0.10 

0.026 

0.180 

0.008 

0.055 

0.01 

0.05 

0.020 

0.137 

0.006 

0.042 

0.05 

0.10 

0.026 

0.183 

0.008 

0.056 

0.01 

0.10 

0.025 

0.172 

0.008 

0.052 

0.01 

0.05 

0.019 

0.131 

0.006 

0.040 


the uncertainties on the Einstein radius and central velocity 
dispersion, namely 

(£is,£o) = (0.05,0.10) , (0.01,0.10) , (0.01,0.05) . 

As a final ingredient, we need to set the total number of 
lenses and the binning. We consider two cases. First, we 
take 1000 lenses splitted in 20 equally spaced bins over the 
redshift range (0.1, 2.1), while secondly we take 10000 lenses 
splitted in 50 bins over the same redshift range as improved 
sample. 


4.2 Results 

Table[T] and [2] summarizes the results of the Fisher matrix 
forecast for different models, sample and error configura¬ 
tions. First, we consider the simplifying assumption that Wa 
is knowrQ. While for ACDM and Quiessence models this 
means that, in a hypothetical fit, we force Wa = 0, in the 
case of Thawing model, one has to se t Wg = —1.58( 1 -|- wo) 
in order to get the correct behaviour (lLindeill2015al i. 

The numbers in Table[T] convincingly show that the dis¬ 
tance ratio method is quite effective at constraining the 
cosmological parameters (Qm,wo) if a prior on Wa is set. 
No matter which model is taken as fiducial, we find that 
1000 lenses are sufficinet to get (t{Q,m)/^m 10% and 

a{wo)/wo ~ 13 — 20%. Increasing the sample by an order 
of magnitude roughly halves these numbers pointing at a 
possible saturation of the method accuracy. For fixed model 
and lens sample, it turns out that reducing the error on go is 
a better strategy to improve the constraints since it is more 
effective at narrow down the error on wo- The constraints 
are basically the same for the three models, although those 
for the Thawing case are consistently smaller. This can be 
qualitatively explained noting that, in this case, wo has a 
larger impact on the distance ratio thanks to the presence 
of the exponential term in the dimensionless Hubble param¬ 
eter (which drops out when Wa = 0 is used). 

Although encouraging, these results are based on the 


^ From the point of view of the Fisher matrix, this mean s we 
delete the row and columns corresponding to Wa dCo 3 l2009ll . 
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Figure 2. 68% contour plots for the fiducial ACDM model as inferred from 1000 (solid) and 10000 (dotted) lenses samples under the 
three different assumptions on the errors (black, blue and red) in the same order as in the text. 


strong prior that we know how to set Wa- In a realistic 
application, however, one only assumes that the DE EoS 
is given by the CPL approximation and fits for the three 
parameters {Q.m,'W o,Wa)- The constraints we get are sum¬ 
marized in Table[2]and show a significant degradation with 
respect to the case with Wa fixed. While the trends with 
{eejEo) are qualitatively the same, the number of lenses 
becomes now of primary importance. Indeed, 10000 rather 
than 1000 lenses are now needed to get a{^M )/^m ~ 10% 
and a{Q.M)/^M ^ 15%, while only weak constraints can 
be put on Wa- This result can be easily understood look¬ 
ing at the pivot redshift Zp defined as the value of z where 
the errors are uncorrelated. Not surprisingly, we find that 
Zp ~ 0.20 no matter which error configuration, lens sample 
or fiducial model is used. This is an obvious consequence 
of the lens redshift distribution which indeed has a median 
value close to the pivot redshift. As a result, the DE EoS is 
best constrained at Zp with the error on Wp = w{zp) turning 
out to b(|3 

a{wp) = (0.081,0.076,0.059) 

from 1000 lenses in the three error configurations and 
(j{wp) = (0.024,0.023,0.018) 

from 10000 lenses. These results clearly show that the dis¬ 
tance ratio method is better suited at probing the low red¬ 
shift behaviour of the DE EoS so that its redshift evolution 
is hardly probed thus explaining the weak constraints on Wa ■ 
It is interesting to look at the FoM of the method. We 
get FoM = (2.0, 2.2, 3.9) when using 1000 lenses, while it 
is FoM = (21.2,23.9,41.9) if 10000 lenses are used. It is 
evident that interesting values can only be obtained reduc¬ 
ing as much as possible the error on (Tq and increasing the 
lens sample. We, however, note that the orientation of the 
{wo,Wa) ellipses (see Fig,[^, is quite different from other 
classical methods such as SNela and CMBR thanks to the 
different pivot redshift. It is therefore likely that a combi¬ 
nation of the distance ratio method with these other probes 
helps breaking degeneracies thus boosting the total FoM. 

® Hereafter, we only give the results for the fiducial ACDM, but 
they are quite similar for other models. 


Table 2. Same as Tablc[T]but for (DjvriuJOirua). 


se 

£0 


o - ifc ( ioo ) 

criki'^a) 

(^lOki^M) 

o-iofc (^ i ; o ) 

(^lOki'^a) 

0.05 

0.10 

0.092 

0.47 

2.58 

0.028 

0.15 

0.79 

0.01 

0.10 

0.087 

0.44 

2.43 

0.026 

0.14 

0.74 

0.01 

0.05 

0.067 

0.33 

1.84 

0.020 

0.10 

0.56 

0.05 

0.10 

0.100 

0.44 

2.48 

0.030 

0.14 

0.76 

0.01 

0.10 

0.095 

0.42 

2.34 

0.029 

0.13 

0.72 

0.01 

0.05 

0.073 

0.31 

1.77 

0.022 

0.10 

0.54 

0.05 

0.10 

0.093 

0.44 

2.38 

0.028 

0.14 

0.73 

0.01 

0.10 

0.088 

0.42 

2.24 

0.026 

0.13 

0.69 

0.01 

0.05 

0.067 

0.31 

1.69 

0.020 

0.10 

0.52 


5 CONCLUSIONS 

Einstein rings have always attracted a lot of attention not 
only for their spectacular beauty, but also as a probe of the 
dark matter content of lens galaxies. The measurement of 
both the Einstein radius and central velocity dispersion al¬ 
lows to strengthen the constraints on the halo density pro¬ 
file, but can also be used as a way to estimate the ratio 
between the distance to the source and that between lens 
and source and hence probe the cosmic expansion. Moti¬ 
vated by the previous literature results, we have proposed 
a novel approach which, although based on the same idea, 
ameliorates the estimate of the distance ratio. On one hand, 
we have adopted a realistic two components model for the 
lens galaxy taking care of both the stellar and DM contri¬ 
butions. This led us to a general formula which permits to 
estimate the distance ratio = Ds/Dia as a function of ob¬ 
servable quantities and a nuisance parameter Te depending 
on the details of the halo model and the lens parameters. 
The analytical formula we derive for Te has allowed us to 
show that its value is approximately constant for lenses be¬ 
longing to the same redshift bin, so that we finally get an 
approximate estimator which helps us to define a likelihood 
function for the estimate of cosmological parameters. It is 
worth stressing that the method proposed does not depend 
on the particular halo model adopted or on the details of the 
and Cvir-M^ir relations adopted. This is a 
further improvement with respect to the standard approach 
which, on the constrary, postulates that all lenses may be 
described by the singular isothermal sphere. 

A well founded and movitated method to constrain cos- 
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mological parameters should be useless if the numbers of 
tracers and/or the requirements on the uncertainties on the 
observable quantities are too demanding. To investigate this, 
we have carried on a Fisher matrix analysis changing both 
the number of lenses and the relative errors on {Re, £o)- It 
turns out that a survey measuring {Re, cro) with (1, 5)% ac¬ 
curacy for 10000 lenses can get FoM ~ 40, while ~ 10% con¬ 
straints on {flM,wo) can be obtained with 1000 lenses only 
if Wa is fixed. It is therefore worth wondering whether these 
requirements are realistic. To this end, we note that present 
day lenses already achieved es 5% and eo 10%, so that 
it is not unrealistic to guess that higher quality images (as 
can be obtained from future large telescopes and satellite 
missions) and improved spectrograph can easily achieve our 
requirements. More demanding is the constraint on the to¬ 
tal number of lenses. As an example, howeve r, one can note 
that the Euclid mission jLaureiis et ^ 2011) is ex pected to 
observe up to 170000 strong lenses ( Collett l2015l l. so that 
asking that both Re and ao are measured for less than 10% 
of them is not a too demanding requirement. 

Although the present day lens samples are far from the 
numbers explored in our Fisher matrix analysis, it could be 
nevertheless interesting to apply the pro posed method to 
real strong lensing data. I Gao et aP (l2015ll have recently as¬ 
sembled a catalog of 118 strong lensing systems and used 
it to constrain the dark energy equation of state under the 
usual assumption of isothermal sphere model. We plan to 
repeat their analysis using our improved approach in order 
to investigate to which extent the constraints on the cos¬ 
mological parameters depend on the assumed lens model. 
Actually, given the small statistics, it is expected that the 
large errors could prevent a conclusive answer, but such a 
study can provide a comparison benchmark for an analysis 
based on simulated samples. 

It is worth noting that the results are always somewhat 
dependent on the lens redshift distribution. Here, we have 
simulated lens samples in such a way that they grossly re¬ 
produce the main features of the SLAGS dataset since this 
is the most used catalog at the moment. However, future 
datasets can have a different redshift distribution, possibly 
changing the method FoM. Although a more detailed anal¬ 
ysis will be carried on in a future work, we can qualitatively 
anticipate that the pivot redshift will likely be close to the 
median survey redshift. It is therefore possible that combin¬ 
ing different samples with different median redshifts could 
help better tracing the DE EoS over a larger redshift range 
thus increasing the FoM. Note that we are here implicitly 
assuming that the median redshift of the full lens sample is 
the same as the median redshift of the subsample with mea¬ 
sured values of {zi, Za, Re, o-q), which can also not be the 
case because of selection effects. On the contrary, since not 
all the lenses contribute to the FoM in the same way, it is 
possible that an optimal strategy can be worked out maxi¬ 
mizing the scientific return of the companion spectroscopic 
survey needed to measure the lenses velocity dispersions. To 
this end, one can follow t he same method ology presented, 
for time delay distances, in iLindeil ll2015bl L 

As a final remark, we stress that, even if suboptimal 
conditions in terms of number of lenses and/or accuracy in 
the measurement of {Re, o-q) are achieved, the distance ratio 
method we have proposed can be efficiently combined with 
other probes helping to break degeneracies and hence boost¬ 


ing the total FoM. Should the present results be confirmed 
by a more detailed analysis of the lens samples properties, 
new light on dark energy will be shed by the strong regime 
of gravitational lensing. 
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